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1.  Introduction 


This  report  summarizes  signal  processing  techniques  used  for  diagnostic  indicators  and  fault 
detection  for  certain  dynamic  components.  Specifically,  the  components  under  consideration  are 
bearings  and  gears.  There  is  a  significant  amount  of  literature  on  the  research  performed  over 
the  last  several  decades  in  the  development  of  diagnostic  techniques  for  these  components.  The 
diagnostic  techniques  can  be  categorized  into  four  fields;  statistical,  time  domain,  frequency 
domain,  and  time-frequency  domain. 

Accelerometers  are  the  primary  sensors  used  in  monitoring  these  components.  Vibrational  data, 
as  a  function  of  time,  is  obtained  by  sampling  the  output  from  these  sensors.  Of  course, 
continuous  sampling/monitoring  of  the  sensor  output  would  be  best.  However,  this  is  not  usually 
practical  due  to  data  storage  limitations;  consequently,  part  of  the  up-front  analysis  must  be  to 
determine  the  minimum/optimal  time  block  size  which  must  be  recorded.  Sampling  rate  is  also 
critical,  where  a  rate  of  ten  times  the  expected  highest  frequency  component  is  considered 
optimal  and  where  the  lower  bound  is  at  least  two  times  the  expected  highest  frequency 
component. 

Due  to  multiple  overlapping  sources  of  vibration,  the  raw  signal  tends  to  “look  like”  random 
noise,  and  determining  the  components’  health  based  on  the  raw  time  measurements  is  generally 
not  possible.  Simple  statistical  techniques  are  used  to  aid  in  the  interpretation  of  the  vibrational 
measurements.  Transformation  of  the  time  measurements  into  the  frequency  domain,  along  with 
a  combination  of  statistical  techniques,  makes  the  assessment  easier.  When  a  fault  is  generated 
in  a  component,  the  vibrational  information  that  is  related  to  the  fault  is  typically  both  nonlinear 
and  non-stationary.  In  most  cases,  techniques  that  incorporate  Time-Frequency  analysis  are 
more  appropriate.  In  order  to  improve  the  probability  of  fault  detection  while  minimizing  the 
false  alarm  rate  in  a  health  assessment,  these  diagnostic  indicators  will  need  to  be  fused  together. 
There  are  many  techniques  that  are  being  explored  to  automate  the  decision  process  and  to 
reduce  the  burden  of  human  interpretation  of  processed  data. 


2.  Bearings 


2.1  Bearing  Faults 

There  are  many  mechanisms  associated  with  bearing  faults.  The  various  types  of  faults  and  the 
ranking  of  the  occurrence  of  some  fault  types  have  been  reported.  Detection  of  these  faults  is  not 
equal,  but  the  difficulty  becomes  more  challenging  as  the  complexity  of  the  bearing’s  application 
and  environment  increases. 
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A  summary  of  analysis  from  the  literature  in  terms  of  deterioration  of  the  bearings  resulted  in  the 
identifieation  of  the  following  modes  of  failure  (1-3): 

•  Fatigue  -  the  degradation  of  the  material  due  to  normal  usage  over  time.  Minute  eraeks 
develop  in  the  bearing  surface  and  eventually  progress  to  the  surface  where  the  material 
will  separate.  Also  known  as  pitting,  spalling  or  flaking. 

•  Wear  -  normal  degradation  caused  by  dirt  and  foreign  particles  causing  abrasion  of  the 
contact  surfaces  over  time  resulting  in  alterations  in  the  raceway  and  ball  bearings. 

•  Plastic  deformation  -  alterations  in  the  contact  surfaces  as  a  result  of  excessive  loading 
while  stationary  or  during  with  small  movements. 

•  Corrosion  -  the  degradation  as  a  result  of  water  or  other  contaminants  in  the  lubrication  of 
the  bearing.  Oxidation  rust  products  are  formed  on  the  surfaces  and  interfere  with  the 
lubrication  and  rolling  operation  of  the  bearing.  The  subsequent  abrasion  results  in  wear, 
flaking  and  spalling. 

•  Brinelling  -  formation  of  regularly  spaced  indentations  distributed  over  the  raceway 
corresponding  to  the  Hertzian  contact  area.  Possible  causes  are  static  overloading  or 
vibration  and  shock  loads  when  in  a  stationary  position.  This  can  lead  to  spalling. 

•  Lubrication  -  the  lack  of  sufficient  lubricant  that  leads  to  skidding,  slip,  increased  friction, 
heat  generation  and  sticking.  This  can  also  anneal  the  bearing  elements  reducing  their 
hardness  and  fatigue  life. 

•  Faulty  installation  -  includes  excessive  preloading  in  either  radial  or  axial  directions, 
misalignment,  tight  fits,  loose  fits  or  damage  in  installation  process. 

•  Excessive  loads  -  self  explanatory. 

•  Overheating  -  self  explanatory. 

•  Seizing  -  self  explanatory. 

In  terms  of  the  frequency  of  occurrence  related  to  these  fault  mechanisms,  table  1  provides  an 
indication  of  bearing  failures  from  the  manufacturing  process,  shipment  and  storage,  and 
installation  and  operation  (2). 
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Table  1.  Summary  of  fault  mechanism  occurrence. 


Cause 

Frequency  of  occurrence  % 

Wear 

51 

Abrasion 

25 

Fatigue 

18 

Corrosion 

8 

Overload 

49 

Deformation 

14 

Break 

12 

Craek 

12 

Flot  run 

11 

Table  2  outlines  the  relative  degree  of  difficulty  in  detecting  the  faults  associated  with  bearings 
in  various  applications  (1). 

Table  2.  Difficult  in  detecting  bearing  faults. 


Machine  Type 

Class  of  Machine 

Degree  of  Difficulty  for  Bearing 
Fault  Detection 

Fans,  Electric  Motors, 

Generators 

Easiest 

1 

Compressors,  Pumps 

Slightly  Complicated 

2 

Industrial  Gearboxes 

Complicated 

4 

Turbines  including  Gas 

Turbine  Engines 

Difficult 

5 

Flelicopter  Transmissions 

More  Difficult 

7 

Specialized  Rotating 

Machinery  with  Extreme  Noise 
Environments 

Most  Difficult 

10 

2,2  Bearing  Defect  Frequency  Equations 

The  spectral  characteristics  of  bearing  faults  have  been  derived  and  defined  in  table  3.  There  are 
five  basic  motions  that  can  be  used  to  describe  the  dynamics  of  a  bearing  movement  and  each 
motion  generates  a  unique  frequency  response.  The  definition  of  these  fault  frequencies  are  as 
follows: 

•  Shaft  rotational  frequency  (fs)  -  the  rotational  frequency  of  the  rotor  or  shaft  is  fundamental 
to  the  movements  of  bearings.  In  a  steady  state  operating  condition,  the  bearing  outer 
raceway  can  be  assumed  to  be  stationary,  while  the  inner  raceway  is  rotating  at  the  speed  of 
the  shaft. 

•  Fundamental  Cage  frequency  (fc)  -  the  rotational  frequency  of  the  cage  of  the  bearing. 

•  Ball  Pass  Outer  Raceway  frequency  (fspo)  -  defined  as  the  rate  at  which  the  balls  pass  a 
point  on  the  track  of  the  outer  raceway. 

•  Ball  Pass  Inner  Raceway  frequency  (fspi)  -  defined  as  the  rate  at  which  the  balls  pass  a 
point  on  the  track  of  the  inner  raceway. 
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Ball  Spin  frequency  (/g)  -  the  rate  of  the  rotation  of  a  ball  about  its  own  axis. 


The  equations  to  calculate  these  frequencies  are  based  on  characteristics  of  the  bearing  as 
outlined  in  table  3.  Figure  1  defines  each  of  the  parameters  {1,  4-6). 


Table  3.  Bearing  detect  frequency  equations. 


Fault  type 

Defect  frequency 

Fundamental  Cage 

II 

■  d  ] 

1 - cosor 

D  J 

Ball  Pass  Outer 

f 

f  -dj-M 

J  BPO  2 

r.  d 

1 - cos« 

D 

Ball  Pass  Inner 

f 

J  BPI  2 

1 - 1 

t/5 

O 

o 

+ 

Ball  Spin 

II 

(d\ 

J 

1-  —  cos  a 

\D )  \ 

where  -  shaft  rotational  frequency 
a  =  contact  angle 

D  -  Pitch  diameter 

d  -  ball  bearing  diameter 
-  number  of  ball  bearings 
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2,3  Signal  Processing  Techniques 

Signal  processing  techniques  are  algorithms  that  are  applied  to  the  signals  obtained  from  the 
sensor  to  gain  further  insight  into  interpretation  of  the  measured  data  set.  Typically  for  dynamic 
components  such  as  bearings  and  gears,  the  data  measurements  are  obtained  from 
accelerometers.  In  general,  the  sensor  output  is  converted  into  discrete  values  as  a  function  of 
time  by  the  data  acquisition  system  that  reads  the  sensor  output.  Analysis  of  the  sensor  data  is 
very  difficult  in  its  raw  form.  Many  signal  processing  techniques  have  been  developed  to 
process  the  data  into  a  form  that  can  be  more  readily  interpreted.  Signal  processing  techniques 
have  been  developed  for  a  broad  range  of  fields  such  as  communications,  imaging,  radar  and 
sonar  applications.  Many  of  these  techniques  are  applicable  to  the  analysis  of  the  vibrational 
data  from  bearings  and  gears.  These  techniques  are  defined  in  the  areas  of  statistical,  time 
domain,  frequency  domain  and  time-frequency  domain  analysis. 


2,3,1  Statistical  Analysis 

Statistical  analysis  is  the  mathematical  science  dealing  with  the  analysis  or  interpretation  of  data. 
Data  analysts  can  use  a  few  straightforward  statistical  techniques  as  means  of  summarizing  the 
collected  data  from  the  sensors.  These  statistical  techniques  are  under  the  area  of  descriptive 
statistics,  a  methodology  to  condense  the  data  in  quantitative  terms. 


Statistical  techniques  that  are  mainly  used  for  alarm  purposes  in  industrial  plants  are  the 
statistical  moments  of  order  two,  three  and  four  (5).  These  are  also  known  as  the  variance, 
skewness,  and  kurtosis.  The  general  equation  for  the  order  of  moment  is  as  follows: 


M 


r 


1 

N 


where  r  is  the  order  of  the  moment 
N  is  number  of  data  value 
i  is  the  index  of  the  data  value 


X  is  the  mean  value  of  the  data  set 


2. 3. 1.1  Mean 


Mean  is  the  most  common  measure  of  a  statistical  distribution.  In  this  case,  mean  is  the 
arithmetic  average  for  a  set  of  measurements. 


-  1 

x  =  —  V  X 

Ntt  ' 


2. 3. 1.2  Variance 

Variance  is  a  measure  of  the  dispersion  of  a  waveform  about  its  mean,  or  is  called  the  second 
moment  of  the  signal  (7). 
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2. 3. 1.3  Skewness 


Skewness  is  the  statistieal  moment  of  the  third  order  normalized  by  the  standard  deviation  to  the 
third  power.  This  moment  indicates  the  asymmetry  of  the  probability  density  function  or  degree 
of  deviation  from  the  symmetry  of  a  distribution  (7,  8). 


/  = 


2. 3. 1.4  Kurtosis 


Kurtosis  is  the  fourth  statistical  moment,  normalized  by  the  standard  deviation  to  the  fourth 
power.  It  is  a  measure  of  whether  the  data  is  peaked  or  flat  relative  to  a  normal  distribution. 
Most  background  signal  measured  by  a  data  acquisition  is  considered  to  have  a  normal 
distribution  in  amplitude.  The  normal  distribution  has  a  value  of  3  (7,  8). 


K  = 


K  = 


iV^(x.-x) 


K  - 


K  = 


2,3,2  Time  Domain  Analysis 

Time  domain  analysis  is  usually  simpler  than  frequency  domain  analysis  and  less  computational. 
Time  domain  analysis  is  based  on  differentiating  good  bearings  from  defective  ones  by 
evaluating  differences  in  statistics  such  as  root  mean  square  (RMS),  peak  value,  crest  factor  and 
kurtosis.  Pattern  recognition  methods  are  used  to  discriminate  between  good  and  bad 
components  (2). 
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2.3.2. 1  RMS 


RMS  is  related  to  the  energy  of  the  signal.  The  presence  of  defects  are  directly  detected  by  the 
increase  in  vibration  level  (7,  8). 


RMS  = 


N 

ZiM 


23.2.2  Maximum  Amplitude  Value 

Maximum  amplitude  value  indicates  the  severity  of  a  bearing  defect  (7). 
2. 3. 2. 3  Peak  Level 

Peak  =  14  (max(x(t))  -  min  (x(t))) 


2. 3. 2.4  Crest  Factor 


Crest  factor  is  a  measure  of  how  much  impacting  is  occurring  in  the  time  waveform.  Impacting 
in  the  time  waveform  may  indicate  rolling  element  wear  or  cavitation  (7). 


CrestFactor  = 


PeakLevel 

RMS 


2,3,3  Frequency  Domain  Analysis 

Frequency  domain  analysis  is  the  scrutiny  of  data  which  has  been  produced  by  transforming  the 
raw  data  from  the  time  domain  into  the  frequency  domain.  In  this  process,  the  time  series  data  is 
transformed  to  sums  of  simpler  trigonometric  functions  such  as  the  sine  and  cosine  functions. 
Through  this  conversion  or  transformation,  the  detection  of  faults  is  generally  easier,  since  it 
improves  the  signal  to  noise  ratio  of  the  data  and  the  fault  signature  tends  to  be  more  obvious. 

Signals  in  general  can  be  divided  into  two  classes:  stationary  and  non-stationary.  Signals  whose 
average  properties  do  not  change  with  time  and  which  are  independent  of  the  particular  sample 
record  used  to  determine  them  are  said  to  be  stationary.  Non-stationary  signals  are  those  whose 
average  properties  change  with  time  {13). 

2 . 3 . 3 . 1  F  ourier  Analysis 

Stationary  signals  are  characterized  by  time-invariant  statistical  properties,  such  as  the  mean 
value  or  autocorrelation  function.  Analysis  of  stationary  signals  has  largely  been  based  on  well- 
known  spectral  techniques  such  as  the  Fourier  transform,  which  identifies  the  constituent 
frequency  components  within  the  signal  {12).  Traditional  data  analysis  methods  such  as  Fourier 
analysis  are  all  based  on  linear  and  stationary  assumptions,  i.e.,  the  signal  to  be  processed  must 
be  linear  and  temporally  stationary  (55).  In  Fourier  analysis,  the  mathematical  process  is  to 
transform  the  time  measurements  into  the  frequency  domain,  where  it  may  be  easier  to  analyze 
and  interpret  the  sensor  measurements. 
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2. 3. 3. 2  Fast  Fourier  Transform 


The  Fast  Fourier  Transform  (FFT)  is  an  effieient  algorithm  for  caleulating  the  diserete  Fourier 
transform.  A  diserete  Fourier  transform  eonverts  a  series  of  time  diserete  measurements  into  its 
eomponents  in  the  frequeney  domain.  There  are  many  forms  for  implementation  of  the  FFT 
algorithms.  The  FFT  provides  the  same  result  as  the  diserete  Fourier  transform,  but  with  the 
main  differenee  being  the  speed  of  the  overall  ealeulation.  For  large  size  of  data  measurements, 
this  ean  reduee  the  eomputational  time  by  several  orders  of  magnitude. 


V-l 

=0 


«=0 

^  =  0,...,A-1 
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2. 3. 3. 3  Cepstrum  Analysis 

Cepstrum  analysis  is  a  tool  for  the  detection  of  periodicity  in  a  frequency  spectrum.  A  number  of 
faults  tend  to  cause  modulation  of  the  vibration  pattern  resulting  from  tooth  meshing,  and  this 
modulation  gives  rise  to  sidebands  in  the  frequency  spectrum.  Only  experience  would  tell  at 
what  point  the  modulation  is  serious  (10). 

The  principal  use  of  the  cepstrum  for  bearing  fault  detection  is  in  detecting  periodicities 
associated  with  bearing  frequency  harmonics  and  associated  sideband  patterns.  The  cepstrum  is 
the  signal  processing  technique  that  takes  the  inverse  FFT  of  the  logarithm  of  the  squared 
magnitude  of  the  Fourier  transform  of  the  measured  signal. 

C,(r)  =  FF7"'(log[S„(/)]) 


2,3,4  Time-Frequency  Analysis 

Time-frequency  analysis  is  an  attempt  to  overcome  some  of  the  shortcomings  of  Fourier 
analysis.  Time-Frequency  analysis  is  a  series  of  signal  processing  techniques  for  analyzing 
signals  which  are  transient  in  nature  or  non-stationary.  The  statistical  properties  of  a  non¬ 
stationary  signal  change  over  time.  The  time-averaging  properties  of  a  non-stationary  signal 
change  over  time,  thus  making  the  time-averaging  approach  adopted  in  the  Fourier  transform 
ineffective.  Time-frequency  analysis  techniques  are  methods  for  analyzing  nonlinear  and  non- 
stationary  data  (31). 

Time  frequency  techniques  decompose  one-dimensional  time-series  signals  into  a  two- 
dimensional  plane  by  exposing  the  time-dependent  variations  of  characteristic  frequencies  within 
the  signal,  thus  presenting  a  valid  and  effective  tool  for  non-stationary  signal  analysis  (12). 

There  are  many  techniques  that  have  been  developed  to  perform  time-frequency  analysis. 
Unfortunately,  these  techniques  generate  artifacts  in  cross-term  products  and  the  user  must  be 
aware  of  how  application  of  the  particular  technique  will  contaminate  the  results.  All  the 
methods  are  designed  to  modify  the  global  representation  of  the  Fourier  analysis,  but  they  all  are 
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limited  in  one  way  or  the  other.  Neeessary  conditions  for  the  basis  to  represent  a  nonlinear  and 
non-stationary  time  series  {31): 

a.  complete  -  guarantees  the  degree  of  precision  of  the  expansion 

b.  orthogonal  -  guarantees  positivity  of  energy  and  avoids  leakage 

c.  local  -  crucial  for  non-stationarity,  for  in  such  data  there  is  no  time  scale 

d.  adaptive  -  adapting  to  the  local  variations  of  the  data  can  the  decomposition  fully  account 
for  the  underlying  physics 

The  following  is  a  list  of  time-frequency  analysis  that  will  be  described: 

a.  Short  Time  Fourier  Transform  (STFT) 

b.  Wavelet 

c.  Cohen 

d.  Wigner-Ville 

e.  Choi-Williams 

f  Zhao-Atlas-Marks 
g.  Hilbert-Huang 
2.3.4.1  STFT 

The  STFT  is  basically  the  windowing  or  dividing  of  the  raw  data  into  frames  and  applying  the 
Fourier  transform.  This  division  can  be  overlapping  or  non-overlapping  of  the  data.  It  is  an 
attempt  to  analyze  the  non-stationary  characteristics  of  the  signal.  The  resulting  two- 
dimensional  signal  is  typically  visually  displayed  as  a  spectrogram  which  represents  the 
magnitude  squared  of  the  STFT  or  power  variation  in  the  signal  over  time. 

STFT  (r,/)  =  J  x{t)g{t  - 


2. 3. 4. 2  Wavelet  Transform 

The  Wavelet  transform  is  a  technique  similar  to  the  STFT.  This  technique  decomposes  the  data 
with  various  functions  that  are  scaled  in  amplitude  and  time.  Wavelet  Transform  is  a  windowing 
technique  with  variable  sized  regions.  The  terms  “dilations  and  translations”  are  the  processes 
applied  to  the  basis  function  or  mother  wavelet  as  it  is  scaled  in  width  and  location  as  applied  to 
the  data  set.  As  with  the  STFT,  it  provides  the  advantage  of  temporal  resolution  in  addition  to 
the  frequency  information.  Wavelet  analysis  allows  use  of  long  time  intervals,  where  more 
precise  low  frequency  information  is  needed  and  shorter  regions  where  high  frequency 
information  is  needed  (28). 
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23.4.3  Cohen 


WT{s,T)  =  -^^x{t)y/*^ — ^  dt 


This  is  a  general  elass  of  proeessing  that  performs  analysis  in  the  time  and  frequeney  domains.  It 
is  an  attempt  to  overeome  some  of  the  problems  assoeiated  with  the  STFT.  It  uses  bilinear 
transformation  to  transform  measurement  data  into  the  frequeney  domain  while  aeeounting  for 
the  non-stationary  aspeet  in  the  measured  data  set.  The  Wigner-Ville,  Choi-Williams  and  Zhao- 
Atlas-Marks  are  speeial  implementations  of  the  Cohen  teehnique. 

00  00 

Q(C/)=  j  J  dt 

-00  -00 

where  A^{r],T^is  the  ambiguity  function 

00  /  \  /  \ 

A^[ri,T^=  ^  x\t +  —  X*  t-—  (it 

-CO  V  ^  /  V  ^  / 

O  (;/,  r)  A  the  kernel  function 

2. 3. 4.4  Wigner-Ville 

Wigner-Ville  distribution  is  a  very  important  quadratie-form  time-frequeney  distribution  with 
optimized  resolution  in  both  time  and  frequeney  domains.  Wigner-Ville  is  eomputed  by 
eorrelating  the  funetion  with  itself,  the  eorrelation  being  a  produet  of  the  funetion  at  a  past  time 
with  the  funetion  at  a  future  time  (52).  The  Wigner-Ville  is  the  Cohen  teehnique  with  the  kernel 
funetion  defined  as  unity. 

0(77, r)  =  l 

W^{t,f)=\x\t  +  —  X*  t —  e‘^"^''dT 

jL  V  '^J  \ 

where  x  is  the  conjugate  of  x 


2. 3. 4. 5  Choi-Willams 


The  Choi-Williams  distribution  uses  the  exponential  as  the  kernel  funetion  in  the  Cohen 
distribution  to  suppress  the  eross-term  produets  in  the  transformation. 


cwyj) 


(  f  . 

X  t  +  —  X*  t —  e  '^’"^^drjdT 

\  1)  \  1) 
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2. 3. 4. 6  Zhao-Atlas-Marks  (Cone-shaped  kernel) 


The  Zhao-Atlas-Marks  distribution  uses  a  time-lag  kernel  as  the  kernel  funetion  in  the  Cohen 
distribution  for  suppression  of  the  eross-term  produet. 


0(t,r) 


sin(^?7^)^-2W 

7rr]T 


ZAM{t,f) 


il 


sin(;z-;7r) 

TrrjT 


( 

* 

t  -\ - 

X 

t  — 

1  2; 

1  2j 

2. 3. 4. 7  Hilbert-Huang  Transform 

The  Hilbert-Huang  Transform  is  another  time-frequeney  analysis  teehnique  that  eombines  two 
proeessing  teehniques:  Empirieal  Mode  Deeomposition  (EMD)  and  the  Hilbert  transform.  The 
EMD  is  an  algorithm  where  the  signal  is  deeomposed  into  a  set  of  funetions  ealled  Intrinsie 
Mode  Eunetions  (IME)  whieh  is  almost  monoeomponent.  The  IME  represent  the  simple 
oseillatory  mode  versus  the  harmonie  output  of  the  Eourier  transform.  EMD  is  empirieal, 
intuitive,  direet  and  adaptive,  with  a  posteriori  defined  basis  derived  from  the  data  (55). 

The  Hilbert-Huang  Transform  is  defined  as  follows  (30): 

2. 3. 4. 7.1  Empirieal  Mode  Deeomposition 

An  IME  is  defined  as  a  funetion  that  satisfies  the  following  requirements: 

•  1 .  In  the  whole  data  set,  the  number  of  extrema  and  the  number  of  zero-erossings  must 
either  be  equal  or  differ  at  most  by  one. 

•  2.  At  any  point,  the  mean  value  of  the  envelope  defined  by  the  loeal  maxima  and  the 
envelope  defined  by  the  loeal  minima  is  zero. 

Therefore,  an  IME  represents  a  simple  oseillatory  mode  as  a  eounterpart  to  the  simple  harmonie 
funetion,  but  it  is  mueh  more  general.  Instead  of  eonstant  amplitude  and  frequeney  in  a  simple 
harmonie  eomponent,  an  IME  ean  have  variable  amplitude  and  frequeney  along  the  time  axis. 
The  proeedure  of  extraeting  an  IME  is  ealled  sifting.  The  sifting  proeess: 

•  1 .  Identify  all  the  loeal  extrema  in  the  test  data. 

•  2.  Conneet  all  the  loeal  maxima  by  a  eubie  spline  line  as  the  upper  envelope. 

•  3.  Repeat  the  proeedure  for  the  loeal  minima  to  produee  the  lower  envelope. 

The  upper  and  lower  envelopes  should  eover  all  the  data  between  them.  Their  mean  is  mi .  The 
differenee  between  the  data  and  mi  is  the  first  eomponent  hi: 

X(t)  -m\  =  h\  . 
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Ideally  the  construction  of  hi  described  above  should  have  made  it  symmetric  and  have  all 
maxima  positive  and  all  minima  negative.  Also,  hi  should  satisfy  the  definition  of  an  IMF.  After 
the  first  round  of  sifting,  the  crest  may  become  a  local  maximum.  New  extrema  generated  in  this 
way  actually  reveal  the  proper  modes  lost  in  the  initial  examination.  In  the  subsequent  sifting 
process,  hi  can  only  be  treated  as  a  proto-IMF.  In  the  next  step,  it  is  treated  as  the  data,  then 

h\-  mu  =  hn  ■ 

After  repeated  sifting  up  to  k  times,  hi  becomes  an  IMF,  that  is 

hi(k- 1)  -  tn\k  =  h\k  . 

Then,  it  is  designated  as  the  first  IMF  component  from  the  data: 

Cl  =  h\k . 


2. 3. 4. 7.2  The  Stoppage  Criteria  of  the  Sifting  Process 

The  stoppage  criterion  determines  the  number  of  sifting  steps  to  produce  an  IMF.  Two  different 
stoppage  criteria  have  been  used  traditionally: 


•  1 .  The  first  criterion  is  proposed  by  Huang.  It  is  similar  to  the  Cauchy  convergence  test, 

and  we  define  a  sum  of  the  difference,  SD,  as 


Then  the  sifting  process  is  stopped  when  SD  is  smaller  than  a  pre-defined  value. 


•  2.  A  second  criterion  is  based  on  the  number  called  the  S-number,  which  is  defined  as  the 

number  of  consecutive  siftings  when  the  numbers  of  zero-crossings  and  extrema  are  equal 
or  at  most  differing  by  one.  An  S-number  is  pre-selected,  and  the  sifting  process  will  stop 
only  if  for  S  consecutive  times  the  number  of  zero-crossings  and  extrema  stay  the  same, 
and  are  equal  or  at  most  differ  by  one. 

Once  a  stoppage  criterion  is  selected,  the  first  IMF,  Ci,  can  be  obtained.  Overall,  Ci  should 
contain  the  finest  scale  or  the  shortest  period  component  of  the  signal.  We  can,  then,  separate  ci 
from  the  rest  of  the  data  hy  X{t)  -c\=  r\.  Since  the  residue,  ri,  still  contains  longer  period 
variations  in  the  data,  it  is  treated  as  the  new  data  and  subjected  to  the  same  sifting  process  as 
described  above. 


This  procedure  can  be  repeated  to  all  the  subsequent  q’s,  and  the  result  is 


^n-\  Cfi  T/i  . 

The  sifting  process  stops  finally  when  the  residue,  Xn,  becomes  a  monotonic  function  from  which 
no  more  IMF  can  be  extracted.  From  the  above  equations,  we  can  induce  that 
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X{t)  =  Y,Ci+rn 

v=l 

Thus,  a  decomposition  of  the  data  into  n-empirical  modes  is  achieved.  The  components  of  the 
EMD  are  usually  physically  meaningful  because  the  characteristic  scales  are  defined  by  the 
physical  data. 

2. 3. 4. 7. 3  Hilbert  Transform 

Having  obtained  the  intrinsic  mode  function  components,  the  instantaneous  frequency  can  be 
computed  using  the  Hilbert  Transform.  After  performing  the  Hilbert  Transform  on  each  IMF 
component,  the  original  data  can  be  expressed  as  the  real  part  in  the  following  form  : 

X{t)  =  Real  =  . 

7=1 


3.  Gears 


Gear  vibration  monitoring  employs  the  same  techniques  as  bearing  vibration  analysis.  The 
development  of  diagnostic  indicators  or  algorithms  is  further  along  for  gears  than  compared  to 
bearings.  As  in  bearings,  it  uses  statistical  signal  and  spectral  analysis  (27). 

Many  of  the  diagnostic  indicators/algorithms  use  spectral  comparison  where  the  baseline  power 
(magnitude  squared)  spectrum  is  taken  under  well  defined,  normal  operating  conditions  with  the 
machine  in  known  good  condition.  This  ‘baseline’  spectrum  is  used  as  a  reference  for 
subsequent  power  spectra  taken  at  regular  intervals  throughout  machine  life  (2d).  Vibrational 
measurements  are  one  of  the  common  sources  of  data  for  gear  evaluation.  This  is  conducted 
through  the  use  of  accelerometers.  Based  on  various  controlled  experiments  the  accelerometer 
location  and  orientation  appear  to  be  critical  in  effectively  detecting  damage  early  (19). 

A  survey  of  world- wide  airworthiness-related  accidents  of  civil  fleet  helicopters  between  1956 
and  1986  was  performed.  Breakdown  of  transmission  related  accidents  by  component  is  outlined 
in  the  table  4.  Gears  have  a  significant  impact  as  the  source  of  major  accidents  as  compared  to 
bearings.  Approximately  19%  of  all  transmission  related  helicopter  accidents  were  caused  by 
gear  failures  (26). 
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Table  4.  Summary  of  helicopter  accidents 


Component 

Percentage  of  accidents 

Tail  rotor  drive  shaft 

31.9 

Gears 

19.1 

Main  rotor  drive  shaft 

14.9 

Lubrication  system 

8.5 

Main  gearbox  input  shaft 

8.5 

Bearings 

4.3 

Freewheels 

4.3 

Cooling  fan  drive 

4.3 

Unknown 

4.3 

3.1  Diagnostic  or  Signal  Enhancement  Techniques 

Interpretation  and  eorrelation  of  vibrational  analysis  results  are  diffieult,  even  for  the  most 
experieneed  personnel.  Many  diagnostie  indieators  or  eondition  indieators  have  been  developed 
over  the  last  few  deeades.  These  indieators  typieally  use  a  eombination  of  statistieal  and 
advaneed  signal  proeessing  teehniques. 

Defeets  or  damage  will  inerease  the  maehinery  vibration  level.  Through  vibrational 
measurements,  the  health  of  the  monitored  maehine  is  eontained  in  this  vibration  signature. 
Henee,  the  new  or  eurrent  vibration  signatures  eould  be  eompared  with  previous  signatures  to 
determine  whether  the  eomponent  is  behaving  normally  or  exhibiting  signs  of  failure.  In 
praetiee,  sueh  eomparisons  are  not  effeetive.  Due  to  the  large  variations,  direet  eomparison  of 
the  signatures  is  diffieult.  Instead,  a  more  useful  teehnique  involves  the  extraetion  of  features 
from  the  vibrational  signature  data.  Ideally,  these  features  are  more  stable  and  better  behaved 
than  the  raw  signature  data  itself.  These  features  also  provide  a  redueed  data  set  for  the 
applieation  of  pattern  reeognition  and  traeking  teehniques. 

3.2  Time  Synchronous  Average  (TSA) 

TSA  is  a  rather  eommon  method  for  early  deteetion  of  failure  in  gears.  By  synehronizing  the 
vibration  signal  with  the  rotation  of  a  partieular  gear  and  evaluating  the  ensemble  average  over 
many  revolutions  with  the  start  of  eaeh  frame  at  the  same  angular  position,  a  signal,  ealled  time 
synehronous  average  is  obtained,  whieh  in  praetiee  eontains  only  the  eomponents  whieh  are 
synehronous  with  the  revolution  of  the  gear.  This  proeess  reduees  the  effeets  of  all  other 
sourees,  ineluding  other  gears,  and  noise. 

TSA  is  a  signal  proeessing  teehnique  that  is  used  to  extraet  repetitive  signals  from  additive  noise. 
In  this  teehnique,  a  signal  that  is  synehronous  with  the  shaft  rotation  provides  the  registration 
mark  that  is  used  to  delineate  the  beginning  data  point  assoeiated  with  every  shaft  rotation.  A 
number  of  data  sets  for  individual  shaft  rotations  are  averaged  together  to  a  single  representation 
of  the  bearing  data  set.  Inereasing  the  number  of  averages  should  improve  the  average  by 
averaging  out  more  of  the  non-synehronous  eomponents  and  enhaneing  the  synehronous 
eomponents,  but  it  is  more  eomputationally  intensive  {16). 
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The  TSA  is  used  in  eonjunetion  with  many  of  the  diagnostic  indicators  as  defined  in  the 
following  text.  Terminologies  used  to  describe  these  algorithms  are  somewhat  vague  and  do  not 
necessarily  have  clear  definitions.  Two  of  these  terms  are  residual  and  difference. 

Residual  consist  of  the  TSA  signal  with  the  primary  meshing,  shaft  components  and  their 
harmonics  removed.  What  is  not  clear  is  how  many  harmonics  to  remove  for  the  primary  mesh 
and  shaft  components.  Favorable  results  have  been  obtained  by  high  pass  filtering  the  data  about 
some  frequency  and  only  removing  the  meshing  frequency  and  all  harmonics.  The  residual  is 
used  in  NA4  and  NA4*  diagnostic  indicator. 

Difference  is  the  removal  of  the  regular  meshing  components  from  the  TSA.  The  regular 
meshing  components  consist  of  the  shaft  frequency  and  its  harmonics,  the  primary  meshing 
frequency  and  its  harmonics,  and  the  first  order  sidebands.  This  is  similar  to  Residual  with  the 
addition  of  the  sidebands  of  the  primary  meshing  frequency.  The  difference  is  used  in  FM4  and 
M6  diagnostic  indicator  {16). 

3.3  FMO 

FMO  is  an  indicator  of  major  faults  in  a  gear  mesh.  Localized  faults  tend  to  increase  in  Peak-to- 
Valley  faster  than  the  magnitude  of  meshing  harmonics.  It  is  defined  as  the  ratio  of  the  peak-to- 
peak  level  of  the  signal  average  to  the  sum  of  the  RMS  levels  of  the  meshing  frequency  and  its 
harmonics  {20,  23). 

The  following  are  some  of  its  characteristics  as  defined  by  R.  M.  Stewart  (the  developer)  (18): 

1 .  It  detects  changes  in  the  average  of  a  significant  kind. 

2.  It  is  a  global  technique  and  will  react  to  changes  occurring  anywhere  in  the  frequency  span 
of  the  average. 

3.  It  will  react  to  a  change  in  only  one  frequency  (a  fairly  common  occurrence  with  gearbox 
faults) 

4.  It  is  non-dimensional  in  a  way  that  makes  it  relatively  insensitive  to  load  changes,  but  not 
speed. 

5.  It  is  also  fairly  tolerant  to  accelerometer  positioning  error  . 

The  equation  for  FMO  is  {16,  17,  20,  22,  24,  26): 
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FMO  = 


where 


[Peak  -to-  Peak) 


TSA 


E4f,) 


i=l 


[Peak  -  to  -  Peak) is  the  peak-to-peak  value  of  the  TSA  waveform 

n 

^  (yi )  is  amplitude  of  the  gear  mesh  fundamental  and  harmonics  in  the 

1=1 

frequency  domain 


3.4  FM4 

FM4  is  an  indicator  of  changes  in  the  vibration  pattern  resulting  from  damage  to  a  single  tooth. 

In  this  routine,  the  regular  meshing  components  are  filtered  from  the  signal  average  obtained 
from  the  TSA  operation  and  two  statistical  operations,  standard  deviation  and  kurtosis  performed 
on  the  difference  signal.  The  difference  signal  is  constructed  by  removing  the  shaft  frequency 
and  its  first  few  harmonics,  the  meshing  frequency  and  its  harmonics,  and  their  first  order  side 
bands  (20,  23). 

An  important  problem  in  gearbox  monitoring  is  how  to  treat  the  gearbox,  which  operates  over  a 
wide  range  of  loads  and  speeds  and  is  not  important  enough  to  warrant  holding  in  store  a 
sufficient  number  of  signatures.  The  concept  behind  the  number  is  that  the  average  of  a  perfect 
gear  will  contain  only  shaft  order  harmonics  that  produce  an  even  pattern  with  the  average,  and 
that  many  of  the  coefficient  values  in  the  FFT  of  the  average  appear  solely  in  order  to  the  make 
one  tooth  different  from  its  neighbors.  If  one,  therefore,  takes  an  average  and  then  subtracts  one 
from  the  ‘regular’  components,  and  then  reconstructs  a  ‘regular’  average,  and  finally  subtracts 
one  from  the  other  (hence  the  name  ‘bootstrap  reconstruction’),  once  might  end  up  with  some 
useful  information.  The  main  problem  is,  of  course,  connected  with  choice  of  harmonics  which 
constitute  the  ‘regular’  set  (18). 

The  equation  for  FM4  is  (16,  17,  19,  20,  22,  24,  26,  27): 

-A 

N'^(d,-d) 

FM4  =  — ^ 

r  ^  /  — n2  1 

E(rf, -d) 

_  i=l 

where 

d  is  difference  signal 
d  is  mean  value  of  difference  signal 
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3.5  NA4 


NA4  is  an  indicator  used  to  detect  the  onset  of  damage,  but  it  also  continues  to  react  to  the 
damage  as  it  increases.  This  is  similar  to  FM4,  but  the  residual  signal  is  constructed  by 
removing  regular  meshing  components  from  the  original  signal  while  keeping  the  first  order 
sidebands  in  the  residual  signal.  A  residual  signal  is  created  by  removing  only  the  meshing 
frequency  components  from  the  vibration  signal  and  dividing  the  fourth  statistical  moment  by  the 
current  run-time-averaged  variance,  i.e.,  the  averaged  variance  of  all  signals  from  the  start  of  the 
run  up  to  the  present  {21,  23). 


The  equation  for  NA4  is  {16,  17,  19,  22,  24,  27): 


NA4  = 


1  M 

Mji 


S(' 

i=\ 


-2 


where 

r  is  residual  signal 

r  is  mean  value  of  residual  signal 

N  is  the  total  number  of  time  points  in  the  time  record 

i  is  data  point  number  in  time  record 

M  is  the  current  time  record  number  in  run  ensemble 

j  is  time  record  index  in  run  ensemble 


3.6  M6 

M6  is  an  indicator  of  fault  in  the  gears.  The  algorithm  is  similar  to  the  application  of  kurtosis  to 
detection  of  gear  faults.  In  this  algorithm,  the  sixth  moment  of  the  difference  signal  is 
normalized  by  the  variance  raised  to  the  third  power  {15,  16,  19,  24). 

M6  = - ^ 

_  ^=1 

where 

d.  is  the  difference  signal 
d  is  the  mean  value  of  difference  signal 
i  is  data  point  index 
N  is  number  of  data  points 
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3.7  NB4 


NB4  is  similar  to  the  NA4  indicator.  NB4  uses  the  same  equation  as  NA4  except  that  the 
residual  signal  is  replaced  by  the  envelope  of  the  signal  bandpassed  about  the  mesh  frequency. 
Using  the  Hilbert  transform,  a  complex  time  signal  is  created  in  which  the  real  part  is  the  band¬ 
pass  signal,  and  the  imaginary  part  is  the  Hilbert  transform  of  the  signal.  The  envelope  is  the 
magnitude  of  this  complex  time  signal,  and  represents  an  estimate  of  the  amplitude  modulation 
present  in  the  signal  due  to  the  sidebands  {21). 


The  equation  forNB4  is  {16,  17,  19,  22,  24,  27): 


NB4  = 


^  /  — n4 

jvX(£,-£) 


1  M 

-E 


i=\ 


-2 


where 

E  is  the  envelope  of  the  band-passed  signal 

m=^m)f+H[A{t)f 

A{t)  is  the  band-passed  signal 

H{A{t))  is  the  Hilbert  trasnsform  of  the  band-passed  signal 


3.8  FM4* 

FM4*  is  an  improvement  on  the  FM4  algorithm.  In  this  case,  the  numerator  is  calculated  the 
same  way  as  in  FM4.  Then,  the  denominator  is  replaced  with  the  run  ensemble  average  as 
defined  in  NA4*  {24). 

3.9  NA4* 

NA4*  is  used  to  overcome  the  associated  problem  of  NA4.  In  gears,  as  the  damage  progresses 
from  localized  to  distributed,  the  variance  of  the  kurtosis  increases.  The  NA4  algorithm  uses  the 
kurtosis.  Since  the  kurtosis  is  normalized  by  the  variance,  this  results  in  the  kurtosis  decreasing 
to  normal  values  even  with  damage.  To  overcome  this  problem,  in  NA4*,  the  data  record  is 
normalized  by  the  squared  variance  for  a  good  gearbox  {24). 

3.10  NB4* 

NB4*  is  an  improvement  on  the  NB4  algorithm.  In  this  case,  the  numerator  is  calculated  the 
same  way  as  in  NB4.  The  denominator  is  replaced  with  the  run  ensemble  average  as  defined  in 
NA4*  {24). 

3.11  NP4 

NP4  is  similar  to  the  other  condition  indicators  in  the  sense  that  it  is  normalized  by  the  variance. 
The  difference  is  that  the  time  domain  power  signal  is  derived  from  the  Wigner-Ville  distribution 
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and  does  not  require  the  eomparison  of  the  undamaged  (baseline)  gear  vibration  signal  with  the 
aequired  gear  vibration  signal.  An  advantage  of  this  technique  is  that  it  avoids  the  need  to 
interpret  the  Wigner-Ville  analysis  (25). 


1  N 

Atp4=— y 

where 
P{t)  =  \\W{tJ)\df 


4.  Summary  of  the  Results  Obtained  by  Researchers  in  the  Application  of  the 
Diagnostic  Indicators 


The  following  results  were  noted  from  the  diagnostic  techniques  in  the  literature: 

•  RMS,  Energy  Ratio,  Energy  Operator,  Kurtosis,  and  NA4  are  very  sensitive  to  torque 
fluctuations  and  thus  may  not  be  effective  {19). 

•  Generally  the  RMS  value  of  the  vibration  signal  is  a  very  good  descriptor  of  the  overall 
condition  of  the  gearbox.  Its  value  increases  as  tooth  failure  progresses.  RMS  does  not 
increase  with  the  isolated  peaks  in  the  signal;  consequently  this  parameter  is  not  sensitive 
to  incipient  tooth  failure  (22). 

•  EM4,  NA4,  NB4,  NA4*,  NB4*  were  found  to  be  more  robust  for  diagnosis  of  initiation  and 
progress  of  damage.  NA4*  is  less  likely  to  ignore  existing  damage  (27). 

•  EM4,  NA4,  NB4,  and  NB4*  are  good  indicators  of  initial  pitting  only  (27). 

•  Most  effective  (in  decreasing  order)  were  M6A*,  EM4*  and  NB4.  They  are  sensitive 
enough  to  pick  up  damage  while  not  being  overly  sensitive  to  torque  fluctuations  {19). 

•  NP4  deteriorates  as  the  severity  of  the  damage  increases  on  multiple  gear  teeth  (25). 

•  Crest  factor  indicates  damage  in  an  early  stage.  As  damage  progresses  the  RMS  value  of 
the  vibration  signal  increases  its  value  and  the  crest  factor  decreases  (22). 

•  Kurtosis  computation  has  an  inherent  characteristic  of  being  self-limited,  which  reduces  the 
scope  of  many  kurtosis-based  applications.  The  key  to  developing  an  effective  application 
lies  in  producing  an  appropriately  tailored  signal  (27). 

•  Wigner-Ville  has  two  problems  related  to  processing  artifacts  such  as  cross-terms  and 
aliasing  (25). 
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•  In  terms  of  the  time-frequeney  analysis,  the  best  resolutions  are  provided  by  Zhao-Atlas- 
Marks,  Choi-Williams,  and  Cohen  in  that  order,  but  require  greater  eomputation  time 
compared  to  the  SIFT  {11). 


5.  Automation  of  Fault  Detection 


The  traditional  methods  for  fault  diagnosis  are  categorized  as  pattern  classification,  knowledge- 
based  inference,  and  numerical  modeling.  Pattern  classification  and  knowledge-based  inference 
techniques  are  used  in  the  industry.  In  these  two  methods,  a  human  expert  looks  for  particular 
patterns  in  the  vibration  signature  that  might  indicate  the  presence  of  a  fault  in  the  bearing. 
Alternatively,  numerical  modeling  which  employs  statistical  analysis  and  Artificial  Neural 
Networks  are  used  in  automated  fault  detection  systems  {14). 

The  majority  of  diagnostic  techniques  have  some  ability  to  indicate  the  presence  of  damage. 

Each  technique  has  associated  strengths  and  limitations.  A  technique  may  be  very  good  at 
detecting  particular  faults  but  weak,  or  useless,  at  detecting  others.  Any  one  technique  is 
generally  insufficient  for  detecting  all  various  types  of  damage.  Model-based  techniques  offer  a 
fundamentally  different  approach  to  detection.  The  fundamental  idea  is  that  the  technique  is 
trained  to  recognize  signals  from  a  healthy  system  and  indicate  when  the  vibration  deviates  from 
this  nominal  condition.  Neural  networks  are  the  most  commonly  used  for  model-based 
diagnostics.  A  neural  network  is  defined  as  a  massively  parallel  distributed  processor  that  has  a 
natural  propensity  for  storing  experiential  knowledge  and  making  it  available  for  use.  It 
possesses  two  fundamental  properties:  (1)  the  network  obtains  knowledge  through  a  learning 
process,  and  (2)  the  intemeuron  connection  strengths  or  weights  are  used  to  retain  the 
knowledge.  In  general,  a  neural  network  damage  detection  system  does  not  process  the  vibration 
signal  directly.  Instead,  it  takes  as  input  the  results  of  various  processing  techniques,  learning 
how  the  techniques  behave  in  the  presence  of  healthy  and  faulty  data  {17). 

Artificial  neural  networks  can  approximate  almost  any  non-linear  function,  which  provides  the 
fundamental  fact  that  it  is  possible  to  model  the  non-linear  dynamics  of  gearboxes  by  an 
appropriate  neural  network  (55). 

Self-organizing  feature  maps  (SOFM)  are  a  paradigm  of  neural  networks  to  map  a  high¬ 
dimensional  space  on  to  a  low-dimensional  feature  map.  This  approach  preserves  with  most  of 
the  topological  relationships  of  the  signal  domain  so  that  humans  can  more  readily  understand 
the  high-dimensionality  data  that  is  typically  hard  to  comprehend.  The  advantage  of  SOFM,  as 
opposed  to  traditional  clustering  methods,  is  to  visualize  the  complicated  cluster  distribution  on 
the  low-dimensional  feature  map  without  setting  the  number  of  clusters  in  advance  (55).  This  is  a 
nonlinear  projection  method  that  efficiently  maps  different  characteristic  features  into  the 
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clusters  on  the  map  without  any  explicit  modeling  of  the  system.  In  this  approach  SOFM  to  be 
effective,  feature  selection  is  critical  (36). 


6.  Summary 


There  has  been  mueh  written  on  the  development  of  diagnostic  techniques  for  dynamic 
components  based  on  decades  worth  of  research.  Most  of  the  measured  data  is  obtained  from 
accelerometers  that  provide  vibrational  information.  The  diagnostic  techniques  can  be 
categorized  into  several  domains:  statistical,  time,  frequency  and  time-frequency  domains. 

Most  of  the  diagnostic  techniques  are  basic  signal  processing  techniques  that  have  been 
developed  for  other  fields  such  as  communication,  radar  and  image  processing.  The  basic  signal 
processing  techniques  are  common  for  both  the  bearing  and  gear  fault  detection.  Development 
of  the  bearing  fault  equations  provides  for  an  easier  fault  detection  process.  Since  the  gear  has 
been  identified  as  the  major  component  that  results  in  accidents,  there  has  been  development  in 
terms  of  specific  signal  enhancements  for  improved  detection  of  gear  faults. 

All  the  researeh  has  been  conducted  only  on  particular  test  situations,  which  do  not  necessary 
cover  all  the  degradation  mechanisms.  To  improve  the  confidence  level  associated  with  the 
detection  of  a  fault,  these  diagnostic  indicators  will  have  to  be  combined.  Artificial  Neural 
Network  is  a  possible  solution  of  these  fault  detection,  since  an  individual  fault  detector  probably 
does  not  have  a  low  false  alarm  probability,  nor  a  sufficient  highly  probability  of  detection  by 
itself. 
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